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We study the decoherence of a renormalised quantum field theoretical system. We consider our 
novel correlator approach to decoherence where entropy is generated by neglecting observationally 
inaccessible correlators. Using out-of-equilibrium field theory techniques at finite temperatures, 
we show that the Gaussian von Neumann entropy for a pure quantum state asymptotes to the 
interacting thermal entropy. The decoherence rate can be well described by the single particle 
decay rate in our model. Connecting to electroweak baryogenesis scenarios, we moreover study the 
effects on the entropy of a changing mass of the system field. Finally, we compare our correlator 
approach to existing approaches to decoherence in the simple quantum mechanical analogue of our 
^ . field theoretical model. The entropy following from the perturbative master equation suffers from 

CO ' physically unacceptable secular growth. 
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O . I. INTRODUCTION 

Recently, we have advocated a new decoherence program [lHl] particularly designed for applications in quantum 
field theory. Similar ideas have been proposed by Giraud and Serreau Q independently. Older work can already be 
interpreted in a similar spirit [tl-fici] . Like in the conventional approach to decoherence we assume the existence of a 
distinct system, environment and observer (see e.g. [HHIEl)- Rather than tracing over the unaccessible environmental 
degrees of freedom of the density matrix to obtain the reduced density matrix /3 re d = Tr^fjo], we use the well known 
idea that loss of information about a system leads to an entropy increase as perceived by the observer. If an observer 
performs a measurement on a quantum system, he or she measures n-point correlators or correlation functions. Note 
that these n-point correlators can also be mixed and contain information about the correlation between the system 
and environment. A "perfect observer" would in principle be able to detect the infinite hierarchy of correlation 
functions up to arbitrary order. In reality, our observer is of course limited by the sensitivity of its measurement 
device. Also, higher order correlation functions become more and more difficult to measure due to their non-local 
character. Therefore, neglecting the information stored in these unaccessible correlators will give rise to an increase 
in entropy. 

In other words, our system and environment evolve unitarily, however to our observer it seems that the system 
evolves into a mixed state with positive entropy as information about the system is dispersed in inaccessible correlation 
functions. The total von Neumann entropy SVn can be subdivided as: 

S vN = S s (t) + S ns (t) = S e s + S E + S ns . (1) 

In unitary theories SVn is conserved. In the equation above S s is the total Gaussian von Neumann entropy, that 
contains information about both the system Sg, environment S E and their correlations at the Gaussian level Sg E 
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(which vanish in this paper), and S ng is the total non-Gaussian von Neumann entropy which consists again of 
contributions from the system, environment and their correlations. Although S v jq is conserved in unitary theories, 
Sg(t) can increase at the expense of other decreasing contributions to the total von Neumann entropy, such as S nK (t). 

In the conventional approach one attempts to solve for the reduced density matrix by making use of a non-unitary 
perturbative "master equation" (l6| . It suffers from several drawbacks. In the conventional approach to decoherence 
it is extremely challenging to solve for the dynamics of the reduced density matrix in a realistic interacting, out- 
of-equilibrium, finite temperature quantum field theoretical setting that moreover captures perturbative corrections 
arising from renormalisation. In fact, we are not aware of any solution to the perturbative master equation that 
meets these basic requirements 1 . Secondly, it is important to note that our approach does not rely on non- unitary 
physics. Although the von Neumann equation for the full density matrix is of course unitary, the perturbative master 
equation is not. From a theoretical point of view it is disturbing that the reduced density matrix should follow from 
a non-unitary equation despite of the fact that the underlying theory is unitary and hence the implications should be 
carefully checked. 



A. Outline 



In this work, we study entropy generation in an interacting, out-of-equilibrium, finite temperature field theory. We 
consider the following action 

S[4>, X] = J d D ^[^, X] = / d^£ o [0] + Coix] + Ant[</>, x) , (2) 

where: 

Co[4>] = -±dp<Kx)d v <t>(x)rr - \m%{t)tf(x) (3a) 

A)M = -\dMd, X {^r -\m 2 x x 2 (x) (3b) 

Ant[<M = -^X 3 (x)- \h X 2 (x)4>{x), (3c) 

where r/^ — diag(— 1, 1, 1, ■ • • ) is the /^-dimensional Minkowski metric. Here, <fi(x) plays the role of the system, 
interacting with an environment xi x )i where we assume that A> /i such that the environment is in thermal equilibrium 
at temperature T. In [l| we studied an environment at temperature T = 0, i.e.: an environment in its vacuum state. 
In the present work, we study finite temperature effects. We assume that (4>) = = (x), which can be realised by 
suitably rcnormalising the tadpoles. 

Let us at this point explicitly state the two main assumptions of our work. Firstly, we assume that the observer can 
only detect Gaussian correlators or two-point functions and consequently neglects the information stored in all higher 
order non-Gaussian correlators (of both <p and of the correlation between (f> and x)- This assumption can of course be 
generalised to incorporate knowledge of e.g. three- or four-point functions in the definition of the entropy Secondly, 
we neglect the backreaction from the system field on the environment field, i.e.: we assume that we can neglect the 
self-mass corrections due to the </>-field on the environment X . This assumption is perturbatively well justified [l[ and 
thus implies that the environment remains in thermal equilibrium at temperature T. Consequently, the countertcrms 
introduced to renormalise the tadpoles do not depend on time too such that we can remove these terms in a consistent 
manner. In fact, the presence of the Ax 3 interaction will introduce perturbative thermal corrections to the tree-level 
thermal state, which we neglect for simplicity in this work. 

The calculation we are about to embark on can be outlined as follows. The first assumption above implies 
that we only use the three Gaussian correlators to calculate the (Gaussian) von Neumann entropy: (4>(x, t)4>(y, t)), 
(■rr(x,t)n(y,t)} and \{{<p{x, t),n(y, t)}). Rather than attempting to solve for the dynamics of these three correlators 
separately, we solve for the statistical propagator from which these three Gaussian correlators can be straightforwardly 
extracted. Starting from the action in equation $Z§, we thus calculate the 2PI (two particle irreducible) effective action 
that captures the perturbative loop corrections to the various propagators of our system field. Most of our attention 
is thus devoted to calculate the self-masses, renormalise the vacuum contributions to the self-masses and deal with 



For example in llTt the decoherence of inflationary primordial fluctuations is studied using the master equation however renormalisation 
is not addressed. In Il8l420ll however, perturbative corrections to a density matrix are calculated in various quantum mechanical cases. 
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the memory integrals as a result of the interaction between the two fields. Once we have the statistical propagator 
at our disposal, our life becomes easier. Various coincidence expressions of the statistical propagator and derivatives 
thereof fix the Gaussian entropy of our system field uniquely [J 0, [l(| . 

In section |TT] we recall how to evaluate the Gaussian von Neumann entropy by making use of the statistical propa- 
gator. We moreover present the main results from In section Hill we evaluate the finite temperature contributions 
to the self-masses. In section ITVl we study the time evolution of the Gaussian von Neumann entropy in quantum me- 
chanical model analogous to equation ^ . This allows to quantitatively compare our results for the entropy evolution 
to existing approaches. In section [V] we study the time evolution of the Gaussian von Neumann entropy in the field 
theoretic case and we discuss our main results. 



B. Applications 

The work presented in this paper is important for electroweak baryogenesis scenarios. The attentive reader will 
have appreciated that we allow for a changing mass of the system field in the Lagrangian (|3a[) : = m^(t). Theories 
invoking new physics at the electroweak scale that try to explain the observed baryon-antibaryon asymmetry in the 
universe are usually collectively referred to as electroweak baryogenesis. During a first order phase transition at the 
electroweak scale, bubbles of the true vacuum emerge and expand in the sea of the false vacuum. Particles thus 
experience a rapid change in their mass as a bubble's wall passes by. Sakharov's conditions are fulfilled during this 
violent process such that a baryon asymmetry is supposed to be generated. The problem is to calculate axial vector 
currents generated by a CP violating advancing phase interface. These currents then feed in hot sphalerons, thus 
biasing baryon production. The axial currents are difficult to calculate, as it requires a controlled calculation of 
non-equilibrium dynamics in a finite temperature plasma, taking a non-adiabatically changing mass parameter into 
account. In this work we do not consider fermions but scalar fields, yet the set-up of our theory features many of the 
properties relevant for electroweak baryogenesis: our interacting scalar field model closely resembles the Yukawa part 
of the standard model Lagrangian, where one scalar field plays the role of the Higgs field and the other generalises to 
a heavy fermion (e.g. a top quark or a chargino of a supersymmetric theory). The entropy is, just as the axial vector 
current, sensitive to quantum coherence. The relevance of scattering processes for electroweak baryogenesis has been 
treated in several papers in the 1990s [2ll - [26| . however no satisfactory solution to the problem has been found so far. 
Quantum coherence also plays a role in models where CP violating particle mixing is invoked to source baryogenesis 
[27M3TI ] . More recently, Herranen, Kainulainen and Rahkila [32|-[34| observed that the constraint equations for fermions 
and scalars admit a third shell at fco = 0. The authors show that this third shell can be used to correctly reproduce the 
Klein paradox both for fermions and bosons in a step potential, and hope that their intrinsically off-shell formulation 
can be used to include interactions in a field theoretical setting for which off-shell physics is essential. The relevance 
of coherence in baryogenesis for a phase transition at the end of inflation has been addressed in (3^413 • 

A second application is of course the study of out-of-equilibrium quantum fields from a theoretical perspective. In 
recent years, out-of-equilibrium dynamics of quantum fields has enjoyed a considerable attention as the calculations 
involved become more and more tractable (for an excellent review we refer to HU). Many calculations have been 
performed in non-equilibrium \<fi 4 (x), see e.g. [39l - l4l| . however also see [i3 - |4^ . The renormalisation of the Kadanoff- 
Baym equations has also received considerable attention j45M47| . Calzetta and Hu [48| prove an TJ-theorem for a 
quantum mechanical 0(A)-model, also see [49[, and refer to "correlation entropy" what we would call "Gaussian 
von Neumann entropy" . A very interesting study has been performed by Garny and Miiller [50l | where renormalised 
Kadanoff-Baym equations in \<j) (x) are numerically integrated by imposing non-Gaussian initial conditions at some 
initial time to- We differ in our approach as we include memory effects before to such that our evolution, like Garny 
and Muller's, is divergence free at to- 

Finally, we can expect that a suitable generalisation of our setup in an expanding Universe can also be applied to 
the decoherence of cosmological perturbations 0, [13, l5lH63| . Undoubtedly the most interesting aspect of inflation 
is that it provides us with a causal mechanism to create the initial inhomogencities of the Universe by means of a 
quantum process that later grow out to become the structure we observe today in the form of galaxies and clusters of 
galaxies. Decoherence should bridge the gap between the intrinsically quantum nature of the initial inhomogeneities 
during inflation and the classical stochastic behaviour as assumed in cosmological perturbation theory. 

II. KADANOFF-BAYM EQUATION FOR THE STATISTICAL PROPAGATOR 

This section not only aims at summarising the main results of [l[ we also rely upon in the present paper, but we 
also extend the analysis of [l[ to incorporate finite temperature effects. 
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There is a connection between the statistical propagator and the Gaussian von Neumann entropy of a system. The 
Gaussian von Neumann entropy per mode Sk of a certain translationary invariant quantum system is uniquely fixed 
by the phase space area A& the state occupies: 



S h {t) 



A fc (t) + 1 



loe 



A fc (t) + 1\ A fc (t) 



■log 



f 



2 ° V 2 J 2 

The phase space area, in turn, is determined from the statistical propagator F ( j > (k,t,t'): 

A 2 k (t) = 4 \F(k, t, t')d t d t ,F{k, t, t') - {d t F{k, t, t')f 



(4) 



(5) 



Throughout the paper we set h = 1 and c = 1. The phase space area indeed corresponds to the phase space area 
of (an appropriate slicing of) the Wigner function, defined as a Wigner transform of the density matrix J2|. For a 
pure state we have Ak = 1, Sk = 0, whereas for a mixed state Ak > 1, Sk > 0. The expression for the Gaussian von 
Neumann entropy is only valid for pure or mixed Gaussian states, and not for a class of pure excited state such as 
cigenstates of the number operator as these states are non-Gaussian. The statistical propagator describes how states 
are populated and is in the Heisenberg picture defined by: 



F$(x;x') 



-Tr 



= 2 Tr 



(6) 



given some initial density matrix operator p(to). In spatially homogeneous backgrounds, we can Fourier transform 
e.g. the statistical propagator as follows: 

F^k,t,t') = J dix-^F^x^e-^ 3 -^ , (7) 

which in the spatially translationary invariant case we consider in this paper only depends on k = \\k\\. Finally, it is 
interesting to note that the phase space area can be related to an effective phase space particle number density per 
mode or the statistical particle number density per mode as: 



n k (t) 



A fc (t)-1 



(8) 



in which case the entropy per mode just reduces to the familiar entropy equation for a collection of n free Bosc- 
particlcs (this is of course an effective description). The three Gaussian correlators are straightforwardly related to 
the statistical propagator: 



(4>(x,t)4>(y,t)) =F (x,t;y,Oli=i' 
(Ti-(x,t)7t(y,t)) = <9 t <9 f <F0(f,i;?7,i')| t=t < 

{{$(x, t), n(y, i)}/2) = dfF^S, t; y, t')\ t=t , . 



(9a) 
(9b) 

(9c) 



In order to deal with the difficulties arising in interacting non-equilibrium quantum field theory, we work in the 
Schwinger-Keldysh formalism (63466| . in which we can define the following propagators: 



iA++(x;x') = Tr[p(ta)T[fo')fo)] 
zA--(.t;.t') = Tr[p(t )T[4>(x')4>(x)] 



iA+-(x;x') 



Tr 



(10a) 
(10b) 
(10c) 
(lOd) 



where to denotes an initial time, T and T denote the anti-time ordering and time ordering operations, respectively. 
We define the various propagators for the field analogously. In equation (|TU|) . «Aj + = iA^ denotes the Fcynman 

or time ordered propagator and *A^~ represents the anti-time ordered propagator. The two Wightman functions 
are given by «A^ + and . In this work, we arc primarily interested in the causal propagator A^ and statistical 

propagator F^, which are defined as follows: 



iA%(x-x') = Tr(fi(to)[4>(x), 
1 



-•Al + (x;aO-*At-(fl 



F^x-x') = -Tr[p(t ){$(x'),4>(x)}] =-(iA-+(x;x') + iA+-(x;x')y 



(11a) 
(lib) 
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We can express all propagators iA^ solely in terms of the causal and statistical propagators: 



iA+-{x;x') = F^x;x')-^%{x;x') 

iA-+{x;x') = F4x;x') + ^iAl(x;x') 

iA++(x;x') = F^x^ + ^sgnit-t'^A^x') 

iA--(x;x') = F^x-x^-hgnit-t'^Alix-x') 



(12a) 
(12b) 
(12c) 
(12d) 



In order to study the effect of perturbative loop corrections on classical expectation values, we consider the 2PI 
effective action, using the Schwingcr-Kcldysh formalism outlined above. Variation of the 2PI effective action with 
respect to the propagators yields the so-called Kadanoff-Baym equations that govern the dynamics of the propagators 
and contain the non-local scalar self-energy corrections or self-mass corrections to the propagators. The Kadanoff- 
Baym equations for the system field read: 

(13a) 
(13b) 
(13c) 



{dl- 


™>l)lAj; + (x]x') - 


f d D y 


iM+ + (x;y)iA++(y;x') 




(x;y)iA~ 4 


- + (y,x') 


= i8 D (x- 


x') 


{dl- 


m l)* A t~( x '> x ') - 




iM++(x;y)tA+-(y;x') 




(x;y)iA~ 4 


-{y;x') 


= 




{dl- 


ml)iA- + {x-x')- ^ 


f d D y 


iM-+(x;y)iA+ + {y-x') 


- iM7~ 

<t> 


(x;y)tA~ 4 


- + {y;x')_ 


= 




{dl- 


m l)* A ^~( x ; x ') - 


f d D y 


iM- + {x-y)iA+-( V] x') 


-iM" 
4> 


(x;y)iA' 4 


-{y;x') 


= -iS D (x 


-x 



where the self-masses at one loop have the form: 

ih 2 , 



iM$ c (x; Xl ) 



iA a x c (x; Xl )y 



iM« c (x; Xl ) = -— (zA~(a;;xi)) -i^ 2 zA~(a;;xi)*A5 (a;;xi). 
Note that we have another set of four equations of motion for the x _n cld. We define a Fourier transform as: 

d^k 



iAf(x;x') = J ^^ t Af(k,t,t'y k <*-^ 
iAf{k,t,t') = J & D -\x - x')iAf(x]x')e- li - {S - s '\ 
such that equations (|13p transform into: 

(d 2 + k 2 + m 2 )tA+ + (k,t,t') + / dh iM++(k,t,t 1 )iA++(k,t 1 ,t') - iM+-{KtM)i^ + {k,h,t') 



(14a) 
(14b) 

(15a) 
(15b) 



(d 2 + k 2 + m 2 ) l A+-(k^t') + 
(d 2 + k 2 + m 2 ) l A^ + (k,t,t') + 



dU 



dti 



iM++ (k, t, h)iA+- (k, t u t') - iM+- (fc, t, h)iA-- [k, h , 
iM7+(k, t, h)iA++ (k, h,t r ) - iMT-(k, t, t')iA7+(k, h , t') 



{d 2 + k 2 +ml)iA--{k 1 t,t')+ / dii iM^ + {k,t,t 1 )tA+-(k,t l ,t')-tM^-(k,t,t 1 )iA^-{k,t 1 ,t') 



(16a) 
-i8{t - t') 
(16b) 

(16c) 

(16d) 
i6(t-t'). 



Note that we have extended the initial time io — > — oo in the equation above. Again, we have an analogous set of 
equations of motion for the x _ field. In order to outline the next simplifying assumption, we need to Fourier transform 
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with respect to the difference of the time variables: 



d D k 



(27T 



iA ab (W) = j d D (x-x')iA c x b (x]x')t 



(17a) 
(17b) 



As already mentioned, we will not solve the dynamical equations for both the system and environment propagators, 
but instead we assume the following hierarchy of couplings: 



(18) 



We thus assume that A is large enough such that the x _n eld is thermalised by its strong self-interaction which allows 
us to approximate the solutions of the dynamical equations for \ as thermal propagators (67j : 



iA++{k») = 



k^k^ + m 2 — le 



+ 27rf(fc M fc" + m£)n?»(|to|) 



iA--(k») = — — : +2^S(k^ + m 2 x )n^(\k°\) 

k^k^ + + le K x 



iA+-(W) = 2Tr5(k^ +m 2 x )[9(-k°) + n^(\k°\)] 
iA-+(k») = 2nS(k^ + m 2 x )[e(k )+n^(\k°\)] , 

where the Bose-Einstein distribution is given by: 

n?(fc°) 



1 



2 /3 fed _ i 



P = 



knT 



(19a) 
(19b) 

(19c) 
(19d) 

(20) 



with ks denoting the Stcfan-Boltzmann constant and T the temperature. Here we use the notation k^k^ = — fc 2 , + k 2 

to distinguish the four- vector length from the spatial three- vector length k = ||«||. We thus neglect the backreaction 
of the system field on the environment field, such that the latter remains in thermal equilibrium at temperature T. 
This assumption is perturbatively well justified [l[. Furthermore, we neglected for simplicity the C(A 2 ) correction 
to the propagators above that slightly changes the equilibrium state of the environment field. Note finally that, in 
our approximation scheme, the dynamics of the system-propagators is effectively influenced only by the 1PI self-mass 
corrections. 

In [![, we have considered an environment field x m its vacuum state at T = and in the present work we 
investigate finite temperature effects. Divergences originate from the vacuum contributions to the self-masses only. 
Since we already discussed renormalisation extensively in [l[, let us just state that the renormalised self-masses are 
given by: 



iM$% n (k, t, t') = (d 2 + k 2 )iZf{k, t, t>) + iM# h (fc, t, , 
where the vacuum contributions Z? h (fc,t,i') to the self-masses are given by: 



Z^(k,t,t') = 



64fc7T 2 

h 2 

64fc7T 2 



7e + log 
7e + log 



k 



2fi 2 \At 
k 



2fi 2 \At\ 



7T 

^ 2 
T ^sgn(At) 



,±4fc|At 



(a(2k\At\) =Rsi(2fc|At|)) 



,±ifcAt 



^ci(2fc| A£|) =f «sgn(Ai)si(2fc|Ai| 



(21) 



(22a) 



, (22b) 



and where iM^ b th (k, t, t') are the thermal contributions to the self-masses that yet need to be evaluated. In deriving 
(|22p . we made the simplifyingassumption m x — > 0. The influence of the environment field on the system field is still 
perturbatively under control [l}. Furthermore, ci(z) and si(z) are the cosine and sine integral functions, defined by: 



ci(z) = - 
si(z) = - 



dt 
dt 



cos(i) 

t 

sin(i) 
t 



(23a) 
(23b) 
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Note that the structure of the self-mass (f2Tj) is such that we can construct relations analogous to equation (|12l) : 



M+-(k,t,t') = MF(k,t,t')--iM$(k,t,t') 

M- + (k,t,t') = M^(k,t,t') + ^iM^(k,t,t r ) 

M++(fc, t, t') = M${k, t, + isgn(i - f)iM$(k, t, t') 

M^~(k, t, f) = Mf(Je, t, f) - \sgn{t - f)iM%(k, t, t') . 



(24a) 
(24b) 
(24c) 
(24d) 



This structure applies to both the vacuum and thermal contributions separately. Thus, Z^ (k,t,t') is the vacuum 
contribution to the statistical self-mass and tZ^(k,t,t') the vacuum contribution to the causal self-mass. Similarly 
we can define the thermal contributions to the statistical and causal self-masses M^ th (k,t,t') and iM^ th (k,t,t'), 
respectively. Of course, we still need to evaluate these expressions. The vacuum contributions follow straightforwardly 
from equation 



z-+(k,t,t') + z+-(k,t,t') 

h 2 r / 

cos(fcAi) 7e + log 



64kw 2 

+ 



2n 2 \At\ 



i(2k\At\) j + sin(fc|Ai|) (si(2A;|At|) - | 



z+-(k,t,t')- z-+(k,t,t') 

h 2 



(25a) 



(25b) 



64fc7T 2 



2 cos(fcAt)sgn(At) ( si(2fc| At\) + - ) + 2sin(fcAt) ( ci(2fe|At|) - te - log 



k 


)] 


_2fi 2 \At\_ 





As before, we are primarily interested in the equations of motion for the causal and statistical propagators, as it 
turns out they yield a closed system of differential equations that can be integrated by providing appropriate initial 
conditions. In order to obtain the equation of motion for the causal propagator, we subtract ()16b[) from (|16c[) and 
use equations pip and (|24| to find: 



(d 2 + k 2 + m 2 )A c 4k,t,t')-(d 2 + k 2 ) J d* 1 Z$(M,<i)AJ(Mi,*')- J t dtiM£ th (M,t 1 )A$(Mi,i / ) = 0. (26) 
In order to get an equation for the statistical propagator, we add equation (|16bj) to (|16cj) . which we simplify to get: 



(d 2 + k 2 + m 2 )F^k, t, f) - (Of + k 2 ) 



dti^(A J t,t 1 )J^(fc,ti,0 - / dhZ^k^t^AKk^t') (27) 

) J — DO 

dtiM& th (M,ti)i^(Mi,t')+ / dtiAf^ th (fc,t,ti)AJ(A ! ,ti J t') = 0. 



Due to the non-locality inherent in any interacting quantum field theory, the "memory kernels" , the memory integrals 
in equation (|27j) above, range from negative past infinity to either t or t' . To make the numerical implementation 
feasible, we insert a finite initial time to by hand and approximate the propagators in the memory kernels from the 
negative past to to with the free propagators inducing an error of the order 0(/i 4 /w4): 



(3 t 2 + fc 2 + m2)F (fc,M') - (df + k 2 ) 



to 



dt 1 Z c 4> (k,t,t 1 )F^ c (k,t 1 ,t')+ dhZ^k^hJF^k^t') (28) 



dhZ$(k,t,h)A 



f ec (k,t llt ')- f dt 1 z£(M > ti)A£(M 1> t / ) 

Jtn 



la 



dt 1 M^ th (fc,t,t 1 )Ff e (k,h,t') - / dtiM£ th (M,*i)i ! *(Mi,t') 



+ I dt 1 M£ th (k,t,t 1 )Af ee (k,t 1 ,t')+ I dt 1 M*] th (k,t,t 1 )A c Jk,h,t>)=0. 
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Here, F^ ee {k,t,t') and A c / ICC (k, t, t') are the free statistical and causal propagators which, depending on the initial 
conditions one imposes at t , should either be evaluated at T = or at some finite temperature. The memory 
kernels need to be included to remove the initial time singularity as discussed in [3, • We postpone imposing initial 
conditions to section [Vj but let us at the moment just evaluate the memory kernels in these two cases. The thermal 
propagators read: 



zA c / rco (fc,M') 



COs(dJ0 (t — t')) 



COsftcU (t-t')) , f l n 



■sin (lj^ (t - t')) 



(29a) 
(29b) 



Here, uj 2 = k 2 + in , where in the case of a changing mass one should use the initial mass. Let us now evaluate the 
"infinite past memory kernels" for the vacuum contributions, i.e.: the memory kernels from negative past infinity to 
to using the two propagators above. The other memory kernels in equation (|28|) can only be evaluated numerically, 
as soon as we have the actual expressions of the thermal contributions to the self-masses. Let us thus evaluate: 



(A? + fc 2 



to 



— oo 
to 



dti^(M > ti)*££(Mi > f)- / dtiZ£(M,f 1 )A2 fa "(M 1> t / ) 



(30) 



32tt 2 



cos[fc(t - tj)] cos^fti - t')} / 1 \ sin[fc(t - tjj\ sin^ftx - t')] 
t — ti uJtj, \2 J t-ti 



64W07T 2 



C0S[ ^i^ *? ] [K + *)(* - *o)] + e ~^ ci IK "*)(*- to)] } 

smh(2PW0) L J 

+ sinMt-tO] r,^ s . + _ + | 

sinh (J/?^) L J 



Due to the fact that the free thermal statistical propagator contains a temperature dependence, the corresponding 
"infinite past memory kernel" is of course also affected. In case we would need the T = vacuum propagators in the 
memory kernels only, one can easily send T — > in the expression above, obtaining the same memory kernels as in 



.1] 



(d? + fc 2 



r di^M.iOifctMi,*')- r dt 1 z£(M,*i)AJ &oe (Mi,0 

J — oo J —oo 



{cos[w0(t - i')]ci [(w + k)(t - t )] + sin[w (i - t')]si + k)(t - t )}} 



(31) 



32W07T- 



III. FINITE TEMPERATURE CONTRIBUTIONS TO THE SELF-MASSES 

In this section, we evaluate all contributions to the self-masses for a finite temperature. 



A. The Causal Self-Mass 

Let us first evaluate the thermal contribution to the causal self- mass. Formally, from equation (|14[) . it reads: 



Ml th (k,At = t-t') = -i M^+ h (k,t,t') - M+- h (k,t,t') 



(32) 



-h 2 



d^fci 



(27T) 



where the superscript F" 1 denotes that we should only keep the thermal contribution to the statistical propagator as 
we have already evaluated the vacuum contribution. The following change of variables is useful: 



d 15 " 1 ^ 1 

(2TT) - 1 ~ (2TT) - 1 J Q 
&D-3 

(27r) D - 1 



dfcifcf 
dfci/cf- 2 



&D-3 



D-2 



dfcxfcf 



{^r- 1 Jo j-i 



dcos(6>)[sin(6>)] 



D-4 



(33) 



2u> 



(2kk 1 ) 



D-3 



[K-^X^-oZ)]^, 
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where we have chosen 9 = Z(fe, fci). In the final line we have changed variables to w = w x (||fc — fci||) = (||fc — k\\ 



m 2 ) 1 / 2 , which clearly depends on 



dimensional sphere S 



D-3. 



Furthermore uj± = (k ± fci)' 



m 2 and denotes the area of the D — 3 



„ -D-2 

2tt— 



r(f-i) 



(34) 



Using equation ([29]) . we have 



Ml th (k,At) = h 2 



Z5-3 



dfciifcf 



dw 



2u> 



(2kk 1 ) D - 
1 



-[K- W 2 )( W 2 -^^ 



{k 2 + m 2 )- 



x — n: 



}2 ) cos \{kl + mJ}5A< sm(wAi) 



(35) 



This contribution cannot contain any new divergences as the latter all stem from the vacuum contribution, which 
allows us to let D 4. Moreover, we are interested, as in section [Til m the limit 
simplifies to: 



0. Equation ([35]) thus 



M£ th (fc,At) 



h 2 sin(fcAi) 

4tt 2 kAt 

h 2 sin(fcAt) 
16^2 kjMfi 



dfci 







27rAt 



sin(2fciAt) 
27rAt 



(36) 



coth 



At coincidence At — > 0, the thermal contribution to the causal self-mass vanishes, as it should. 



B. The Statistical Self-Mass 



The thermal contribution to the statistical self-mass is somewhat harder to obtain. It is given by: 
ML h (k, At = t-t>) = i \M7+(k, At) + M+- h (k, At) 



2 [ <£,th* 
~ 2 J (27T) - 1 



F x (k u At)F x (\\k - At) - i A^(fc 1 , At)A x (\\k - h\\,At) 



(37) 



th 



where of course we are only interested in keeping the thermal contributions. The second term in the integral consists 
of two causal propagators that does not contribute at all at finite temperature. It is convenient to split the thermal 
contributions to the statistical self-mass as: 



M F th (fc, At) = Af£ vac _ th (fc, At) + A/£ th _ th (fc, At) 



where, formally, we have: 



Mi 



,vac— th 



(k,At) = -- 



Af^ th _ th (ft,At) 



tf_ r d - 1 ^ cos(^ x (fc 1 )At)cos(^ x (Hfc-fc 1 |l)Af) 
J J (2n)°-i w x (Ai) Wxdlft-jfciH) : 

h 2 f d - 1 ^ cos(w x (£;i)At) 1 cos(cJ x (||fc - A?i||)Ai) 



1 

3l 8w x (fci) _ i 



1 



(38) 

(39a) 
(39b) 



Here, M^ vac _ th is the vacuum-thermal contribution to the statistical self-mass and M^ th _ th is the thermal-thermal 



contribution. As before, we let D 
Equation (|39ap thus simplifies to: 



M. 



-th 



(A, At) 



{2tt) d - 1 w x {ki) e ^(fei)-i Wx (||fc- jfciH) e^x(llfc-fall) -1 

nass and A^ tll _ th is the thermal-i 
0. Let us firstly evaluate the vacuum-thermal contribution. 

cos(fciAt) . . 



4 and 



87T 2 fcAt J 



dki- 



[sin((fc + h)At) - sinflfc - fci|Ai)] 



(40) 
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We have to take the absolute values in the equation above correctly into account by making use of Heaviside step- 
functions and we can moreover expand the exponential to find: 



M. 



</>,vac— th 



(fc,At) 



!6n 2 kAt 



/ dfci V e~ pnkl I sin ((2fci + k)At) + 29{k 1 - k) sin (kAt) 
J o n=i 1 



(41) 



+ {0(k - fci) - 6>Oi - fc)} sin ((2fci - k)At) 



Integrating over k\ and collecting the terms we get 

h 2 



Ml vac . th (k,At) 



8ir 2 kAt 



J2{ cos(/cAt)(l 



n=l 



-fink 



2At 



{/3n) 2 + (2Aty 



-+sin(fcAi)c 



-Pnk 



pn 



The sum can be performed, resulting in 

h 2 



Ml vac - th (k,At): 



l6ir 2 k(At) 2 



±-. 



sin(fcAt) I — log(l-e-^ fe ) +e"^][] — 



1 2iA t - 



/3n (^n) 2 + (2At) 2 



2iAt 2iAi 



(42) 



•cos(fcAt)< - 



1 / 2?rAt 



P 



coth 



2?rAt 



-1 -e 



-0k 



1 lAf 



^ / 2iAt „ , 2zAt 

^life'' 1 M± T i2± T i8 ' 



P 



, (43) 



where 2-F1 is the Gauss' hypergeometric function. For convenience we quote the low temperature ((3k 3> 1) and the 
high temperature (/3k <C 1) limits of this expression. In the low temperature limit (|43[) reduces to: 



^ vac _ th (fc,At) 



/3fe»l 



cos(fcAt) 



167T 2 fc(At) 2 



2?rAt 



P 



coth 



27rAt 



1 



(44) 



+ c 



-0k 



and its coincidence limit is finite: 



lim Ml 

At->0 



(i),vac— th 



cos(fcAt) TTliw + sin( ^ ) 

/i 2 (tt 2 - 6e"^) 



/2At 



2 At IP 



\ 13 l + (2At/!3)< 



247r 2 fc/3 2 



In the high temperature limit (|4"5)l reduces to: 



^vac-th(A.Ai) 



/3fc<l /l 2 



4tt 2 ^ 



cos(fcAt) J log(^fc) + 7s -l + i^^(l± 



sin(fcAt) 
4fcAt 



2iAt 



(45) 



(46) 



2iAA 



2j E 



There is a mild logarithmic divergence, M^ vllc _ th oc log(f3k), in the limit when /3k —> 0. Also note that when we 
derived equation (|46[) above, we tacitly assumed that also At/ P <C 1. We however only use equation (|46|) to calculate 
the coincidence limit At — » of the statistical self-mass in which case this approximation is well justified: 



lim Ml 

At->0 



vac — th (k, At) 



/3fc<i h 2 



4tt 2 /3 



(log(/3fc) - 1) . 



(47) 



The final remaining contribution to the statistical self- mass is M^ th _ th (k, At) in equation (|39b|) and is much harder 
to obtain. In fact, it turns out we can only evaluate its high (/3k <C 1) and low ((3k 3> 1) temperature contributions 
in closed form. For that reason we present the calculation in appendix [5] and in the current section only state the 
main results. The low temperature (f3k 3> 1) limit of M^ th _ th (k, At) is given by: 



167T 2 fc 



cos(kAt) 



2 7 rAtcoth( 2 ^ 



0k 



2 (50 2 + ll(At) 2 ) 



sin(fcAt) 



27TCOth ( 



/3(/3 2 + (At) 2 ) /3 2 + (At) 2 (0 2 + (At) 2 ) 2 (0 2 + (2At) 2 ) 
fcAt 2/3At 8At 



/3 2 + (At) 2 (3 2 + (At) 2 (/3 2 + (At) 2 ) 2 0(/3 2 + (2At) 2 ) 



(48) 
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Note that this expression is finite in the limit when At — ^ 0: 



lim M£ th _ th (fc,At) - 



At-yQ 



, 3 + /C/3 
87T 2 /fc/3 2 ' 



(49) 



The high temperature (/3k <C 1) limit yields: 



0,th-th 



(A, At) 



167T 2 fc/3 2 

/3sin(£;At) 
+ At 

-k/3±ikAt 



k 2 . i „ „ ,„ (fcAt) 2 d / (fcAt) 2 

y-4 ( 7 E-ci(|fcAt|) + log(|fcAt|)) + ^-^— 2 F 3 (l, 1; 2, 2, l+ 7 ; 



ci(2|fcAt|)-7E-log(^ 
/ 2iAt 2iAt hp 



Clearly, the limit At — > of the self-mass above is finite too 

fe/3<i h 2 



/3cos(fcAt) / ,. „ . .. 7r 



(fc/3) 2 / 2iAt 2iAt up 

^-— J - 2 Fi 4,2 T ;3T ;e"~ 

12 V P P 



(50) 



lim„A/fth-th(fc,At) 



At->0 



327T 2 fc^ 2 



7T 2 + 4fc/3 log 



(51) 



where we ignored the sublcading term in equation (|50|) to derive the coincidence limit above. 



IV. ENTROPY GENERATION IN QUANTUM MECHANICS 

Now the stage is set to study entropy generation in our quantum field theoretical model, let us digress somewhat and 
study entropy generation in the analogous quantum mechanical model first. This allows us to quantitatively compare 
the evolution of the entropy resulting from the perturbative master equation and in our correlator approach. A 
comparison in field theory is not possible so far, due to the shortcomings of the conventional approach to decoherence 
using the master equation as discussed in the introduction. Let us consider the quantum mechanical system of N + 1 
simple harmonic oscillators x and q ni 1 < n < N, coupled by an interaction term of the form h n xqz: 



L E + L 



1 



SE 



2 2\ 
0J Q X ) 



A i i 

1 ( -2 , ,2„2\ 1 u „,„2 



(tin - - 2 hnXq « ' 



(52) 



which indeed is the quantum mechanical D = 1 dimensional analogue of the Lagrangian density in equation ([3]) 
considered before. Here, uj and {w ra } are the frequencies of the oscillators as usual. The x oscillator is the system 
in a thermal environment of {q n } oscillators. We absorb the mass in the time in our action, and the remaining 
dimcnsionlcss mass dependence in the {q n }- 



A. The Kadanoff-Baym Equations in Quantum Mechanics 



The free thermal statistical and causal propagator in quantum mechanics read: 

,. cos(w„ (t — t')) , , . 

A" (t,t') = — sin(wn(t - t 7 )) ■ 



(53a) 
(53b) 



The statistical and causal self-energies of the a;-system at lowest order in perturbation theory are defined by: 



<(M') = -^^[K+-(t,t')) 2 + (*A-+(t,t')) ; 



N 



ihl 



n=l 



M c x (t,t<) = - J2^r ( tA t"(*.* / )) -(*A- + (t,o) J 



(54a) 
(54b) 
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and are calculated as: 



N h 2 

n=l 
JV 

■E 

n=l 
AT 



2 



(F to (t,0) a -4(A^(*,0) 2 



16o;2 



coth 



+ 1^ cos (2w„(i - t')) + coth 2 ( 



N 



M*(i,t') 



]T h 2 n F qn (i, t')A c qn (t, t') = J2^ sin ( 2 ^(t - f )) coth 



n=l 



2 



- 1 



(55a) 
(55b) 



As in our field theoretical model we neglect the backreaction from the system on the environment. The Kadanoff-Baym 
equations for the cc-system for the statistical and causal propagators are now given by: 



{d 2 t +uil)F x (t,t>) 



<ttiM s f(t,t 1 )A°(* 1> t / )- / dhM^t^F^tut') 



(56a) 
(56b) 



where {t, t'} > to = 0. An important difference compared to the held theoretical Kadanoff-Baym equations in (|26|) 
and (|28| is that we do not have to renormalise them. Also, we do not have to consider any memory effects before 
to = 0. We can now straightforwardly solve the Kadanoff-Baym equations above by numerical methods to find the 
statistical propagator and hence the quantum mechanical analogue of the phase space area (O and entropy @ as 
functions of time. 



B. The Master Equation in Quantum Mechanics 

In order to derive the perturbative master equation for our model, we follow Paz and Zurek [l6j . The perturbative 
master equation is obtained straightforwardly from the Dyson series, truncated at second order, as a solution to the 
von Neumann equation and reads: 

1 1 N h 

/3 red (i) = -[Hs(t),Mt)] + -^^[{£(t))St t Mt)] (57) 

n=l 

N , t 

-V / dixK^^tO^Jx^-^^ped^+^W^tx)^,^^-*),^^)}], (58) 
m— 1 

where we follow the notation of Paz and Zurek and define the coefficients: 

i^i(Mi) = ^<{£(t),£(ti)})-^<m<&(*i)> (59a) 

Kit(t,h) = M«<[£(f),&( tl )]>. (59 b) 

Also, note that x(t) = x cos (wo £)+£>.,. sin(woi) / 'uq due to changing back to the Schrodingcr picture from the interaction 
picture (l6j |. The master equation above reduces further to: 

i . r* 

Pred(t) = -[H S (t),p Ied (t)} - / dt!V(ti)[x, &ed(t)]] - ir](t 1 ) [x, {x(~t 1 ) , p rcd (t) }] . (60) 

1 Jo 

In the equation above, we dropped the linear term in equation (|57[) as a time dependent linear term will not affect the 
entropy Q- The noise and dissipation kernels v{t) and r](t) are straightforwardly related to Knm{t, t\) and K^m(t, t\), 
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respectively, and read at the lowest order in perturbation theory [18|, |68[ : 

TV 

1/s. .9 



"(*) = E 



ra=l 



I6u;? 



( coth 2 



+ l cos (2w„t) + coth 



-Mf(t,0) 



»?(*) = E g^- sin (2w„i) coth 



n=l 



(6la) 



(6lb) 



Note that we can easily relate the noise and dissipation kernels that appear in the master equation to the self-mass 
corrections in the Kadanoff-Baym equations. This is an important identity and we will return to it shortly. One thus 
finds: 

/WO = -i[H S {t) + ^n 2 {t)x 2 ,p Icd (t)} - z 7 (t)[x,{p ;c ,p rcd (t)}] - D(t)[x, [x,p ied (t)}) - f(t)[x, [p x ,Pred(i)]], (62) 

where the frequency "renormalisation" £l(t), the damping coefficient j(t) and the two diffusion coefficients D(t) and 
f{t) are given by: 



Qr{t) = -2 / dii77(ii)cos(cj ii) 

JY 



(63a) 



E 



■ coth 



7(t) 
D(t) 



d*iT7(ti) 



sin(wo^ 



{— 2w n (l — eos(oJot) cos(2w„t)) + wo sin(woi) sin(2w n i)} 
/3w„\ f sin([o;o - 2w„]i) sin([w + 2o;„]t) 



cjq — 2w„ 



cj + 2w ra 



dii^(ti) cos(o;oti) 



(63b) 
(63c) 



N 

E 



l6o;2 



coth 



f3to ri 



I 



+ coth 2 



sin([o;o — 2u n ]t) sin([ojo + 2u! n ]t) 
ujq - 2uj n uj + 2w„ 

s'm(ujQt) 

UJ 



/(*) = - / dtii/fa) 
Jo 



sin(a;oii) 



A' 



= E 



-hi 



16a;, 2 , w 



wo 



coth 



(63d) 



{o;o(l — cos(woi) cos(2w„i)) — 2ui n sin(woi) sm(2w„i)} 



loI - 4w 2 



+ coth 2 



/3u n \ _ \ 1 - cos(u;ot) 
2 / y cj 



We are now ready to solve the master equation (|6"2"|) . As we are interested in the evolution of 2-point functions, let us 
make a Gaussian ansatz and project this operator equation on the position bras and kets as follows: 



p Ic d(x,y;t) = (x\p le d(t)\y) =A/ r (t)exp [-a(t)x 2 -a*(t)y 2 + 2c(t)xy] 



(64) 



It turns out to be advantageous to directly compute the time evolution of our three non-trivial Gaussian correlators. 
Analogous methods have been used in [691 ] to analyse decoherence in an upside down simple harmonic oscillator. We 
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can thus derive the following set of differential equations [|| : 



d(P) ^ -2(±{*,*}) (65a) 



dt 4(5 R - c) 2 \ 2 

d(P 2 > ,,..2,^/1 



dt v u ' \2 



2{lo* + fi 2 ) ( -{x,p} ) - 4 7 (t)(^) + 2D{t) (65b) 



dt v u /v ' ^ ' J v ' ,w \2 

These equations are completely equivalent to the master equation when one is interested in the Gaussian correlators 
only. Initially, we impose that the system is in a pure state: 

(i 2 (io)> = ^- (66a) 
(p 2 (io)> = y (66b) 
i{i,p}(t )^ = 0. (66c) 

We can then straightforwardly find the quantum mechanical analogue of the phase space area in equation ([5]) and the 
von Neumann entropy for the system in equation (U). 

Let us finally remark that our Kadanoff-Baym equations (|56p can also be obtained starting from the Feynman- 
Vernon Gaussian path integral exponential obtained in fl8l ] that is normally used to derive the perturbative master 
equation. For example, after integrating out the environment at one loop order (which is usually a first step in deriving 
a master equation), the 1PI equations of motion can be obtained from the effective action: 

/>oo />oo />ti 

S s [x + ]- Ss[x~]- dhTlx+ih)) -x'ih)} + i / dt : / dt 2 [x + (ti)-x-(t 1 )}v{t 1 -t2)[x + (t 2 )-x-(t2)] 
Jo Jo Jo 

3 ftl 

Ah &t2[x + (t l )-x-(t 1 )]i 1 (t l -t2)[x + (t2)+x-(t2)] (67) 
Jo 

/> oo -i />oo />oo 

= Ss[x + }~S s [x-}~ dhTlx+it^-x-ih)}-- V ab / dtx / dt 2 x a (t 1 )iM ab (t 1 ;t 2 )x b (t 2 ). 
Jo z a,b=± Jo Jo 

Here, M ab are the self- masses that can be read off from equations (f55j) and (|24|) . S's[a; ± ] is the free action defined by 
equation (|52[) and r\ and i>, or, equivalently, the causal and statistical self- masses, are given in equation (|61[) . In the 
equation above, T denotes the tadpole contribution, which does not affect the entropy [2|, and reads: 

n— 1 x ' 

which is easily inferred from the interaction term in (|52[) and (|53al) . The quantum corrected equation of motion for 
x(t) follows straightforwardly by variation of equation ([67)) with respect to x ± (t), and setting a; ± (t) equal to x(t). 
More generally, if one would introduce non-local sources for two-point functions in the Feynman- Vernon path integral, 
one would obtain the Kadanoff-Baym equations in ([5! 



C. Time Evolution of the Entropy in Quantum Mechanics 



Lot us discuss our results. It is important to distinguish between the so-called resonant regime and non-resonant 
regime Q. In the former, we have that one or more ui n ~ wo/2, with 1 < n < N. In the non-resonant regime 
all environmental frequencies differ significantly from wo/2 and are as a consequence effectively decoupled from the 
system oscillator. If one wants to study the efficient decoherence of such a system, the non-resonant regime is not the 
relevant regime to consider. 

In figure [1] we show the Gaussian von Neumann entropy resulting from the Kadanoff-Baym equations and from 
the perturbative master equation as a function of time in black and gray, respectively. At the moment, we consider 
just one environmental oscillator N = 1 in the non-resonant regime. Here, the two entropies agree nicely up to the 
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Figure 1: Entropy as a function of time for N = 1 in 
the non-resonant regime. The Gaussian von Neumann 
entropy (black) agrees with the entropy from the mas- 
ter equation (gray) up to the expected perturbative cor- 
rections. The dashed line indicates full thermalisation of 
x. We use uji/uja = 2, h/uj'o = 1 and f3ujo = 1/2. 
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Figure 2: Entropy as a function of time for N = 1 in 
the resonant regime. The Gaussian von Neumann en- 
tropy (black) shows a stable behaviour, unlike the entropy 
from the perturbative master equation (gray) that reveals 
unphysical secular growth. The dashed line indicates full 
thermalisation of x. We use u>i /loo = 0.53, h/uil = 0.1 
and fiojo = 1/2. 
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Figure 3: Entropy as a function of time for N = 50 in 
the non-resonant regime. The Gaussian von Neumann 
entropy (black) agrees with the entropy from the mas- 
ter equation (gray) up to the expected perturbative cor- 
rections. We use uj n /u) e [2, 4], 1 < n < N, h/uj'i = 1/2 
and j3ujo = 1/2. 
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Figure J^: Entropy as a function of time for N = 50 in 
the resonant regime. The Gaussian von Neumann en- 
tropy (black) yields a stable behaviour in time, unlike 
the entropy from the perturbative master equation (gray) 
that reveals secular growth. We use u n /ujo G [0.5,0.6], 
1 < n < N, h/uo = 0.015 and [3u) Q = 1/2. 



expected perturbative corrections due to the inappropriate resummation scheme of the perturbative master equation 
to which we will return shortly. However, let us now consider figure [2] where we study the resonant regime for N = 1. 
Clearly, the entropy resulting from the master equation breaks down and suffers from physically unacceptable secular 
growth. The behaviour of the Gaussian von Neumann entropy from the Kadanoff-Baym equations is perfectly stable. 
Moreover, given the weak coupling h/cuy = 0.1, we do not observe perfect thermalisation (indicated by the dashed 
black line). 

If wc consider ./V — 50 environmental oscillators, the qualitative picture does not change. In figure [3] wc show the 
evolution of the two entropies in the non- resonant regime, and in figure |4] in the resonant regime. The entropy from 
the perturbative master equation blows up as before, whereas the Gaussian von Neumann entropy is stable. In figure 
[3] we randomly select 50 frequencies in the interval [2,4] which is what we denote by w n /wo <E [2,4]. In the resonant 
regime we use w„/wo €[0.5, 0.6]. The breakdown of the perturbative master equation in this regime is generic. 

Just as discussed in [3|, energy is conserved in our model such that the Poincare recurrence theorem applies. This 
theorem states that our system will after a sufficiently long time return to a state arbitrary close to its initial state. 
The Poincare recurrence time is the amount of time this takes. Compared to the N = 1 case we previously considered, 



16 



we observe for N = 50 in figurc[4]that the Poincare's recurrence time has increased. Thus, by including more and more 
oscillators, decoherence becomes rapidly more irreversible, as one would expect. If we extend this discussion to field 
theory, where several modes couple due to the loop integrals (hence N — > oo), we conclude that clearly our Poincarc 
recurrence time becomes infinite. Hence, the entropy increase has become irreversible for all practical purposes and 
our system has (irreversibly) decohered. 

In decoherence studies, one is usually interested in extracting two quantitative results: the decoherence rate and 
the total amount of decoherence. As emphasised before, we take the point of view that the Gaussian von Neumann 
entropy should be used as the quantitative measure of decoherence, as it is an invariant measure of the phase space 
occupied by a state. Hence, the rate of change of the phase space area (or entropy) is the decoherence rate and the 
total amount of decoherence is the total (average) amount of entropy that is generated at late times. This is to be 
contrasted with most of the literature |l5( where non-invariant measures of decoherence are used. The statement 
regarding the decoherence rate we would like to make here, however, is that our Gaussian von Neumann entropy and 
the entropy resulting from the master equation would give the same result as their early times evolution coincides. 
The master equation does however not predict the total amount of decoherence accurately. In the resonant regime the 
entropy following from the perturbative master equation blows up at late times and, consequently, fails to accurately 
predict the total amount of decoherence that has taken place. Our correlator approach to decoherence does not suffer 
from this fatal shortcoming. 



D. Deriving the Master Equation from the Kadanoff-Baym Equations 

The secular growth is caused by the perturbative approximations used in deriving the master equation (|65p . The 
coefficients appearing in the master equation diverge when ui n = ljq/2 which can be appreciated from equation (|63l) . 
However, there is nothing non-pcrturbative about the resonant regime. Our interaction coefficient h is still very small 
such that the self-mass corrections to ujq are tiny. 

Here we outline the perturbative approximations that cause the master equation to fail. In order to do this, we 
simply derive the master equation from the Kadanoff-Baym equations by making the appropriate approximations. 
Of course, equation (|65a[) is trivial to prove. The Kadanoff-Baym equations are given in equation ((56]) and contain 
memory integrals over the causal and statistical propagators. We make the approximation to use the free equation of 
motion for the causal propagator appearing in the memory integrals according to which: 

(S t 2 +w 2 )A^( M ') = 0. (69) 
This equation is trivially solved in terms of sines and cosines. Let us thus impose initial conditions at t = t' as follows: 

A c x {t,t')~ A x ^ cc (t,t') = cos(cj At)A x (t' ,t r ) + —sm(uj At)dtA x (t,t%^ (70) 

Wo 

= sin(woAi) . 

Wo 

Here, At = t — t'. We relied upon some basic properties of the causal propagator (see e.g. equation ([75)1 in the next 
section). Likewise, we approximate the statistical propagator appearing in the memory integrals as: 

F x {t,t') ~ F* Iee (t,t') = eos(oj At)F x (t',t') + — sm(u>oAt)d t F x (t,t%=t> (71) 

w 

= cos(w At)(x 2 (t')> + — sin(w At) ( -{x,p}(t') 
w \2 

where we inserted how our statistical propagator can be related to our three Gaussian correlators, from the quantum 
mechanical version of equation ((9]). Note that expression (|71[) is not symmetric under exchange of t and t' , whereas 
the statistical propagator as obtained from e.g. the Kadanoff-Baym equations of course respects this symmetry. 

Now, we send t' t in the Kadanoff-Baym equations and carefully relate the statistical propagator and derivatives 
thereof to quantum mechanical expectation values. From equation (|56a|) . where we change variables to r = t — t\, it 
thus follows that : 

d?F x (t,t')\ t=t , = -oj 2 (x 2 (t)) - f drM^(r,0) Sm{uJoT) +(x 2 (t)) f drM c x (r, 0) cos(w r) (72) 

Jo w Jo 



l{x,P}(t)\J o drA4 c (r,0)- 



W() 
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Using equations (|61j) and (|63|) . equation ([72)1 above reduces to (|65c[) 



df 



>o + ^ 2 )(i 2 > + (P 2 > - /(*) - 27(t) • 



Here, we used the identities derived in equation (|6ip that relate the noise and dissipation kernels of the master equation 
to our causal and statistical self-mass. In order to derive the final master equation for the correlator (p 2 ), we have to 
use the following subtle argument: 

d?d t ,F(t,t% =t , = ~(f(t)). (73) 

In order to derive its corresponding differential equation, we thus have to act with df on equation (|56a|> and then 
send t' — > t. As an intermediate step, we can present: 

d : (f(t)) = -u%(±{x,p}(t))- / i drA/f(r,0)cosKr) (74) 



2 dt u \ 2 

+ / drM c x (r,0) 



-u sm(cj T)(x 2 (t')) +co S (lo t) (±{x,p}tf)) - dt , {d t F x (t,t% =t ,} 

\ 2 / cj 



where we still have to send t' — > t on the second line. Now, one can use: 

d v {d t F x (t,t')\ t=t ,} = (p 2 {t'))-u 2 (x 2 {t')) . (75) 
In the light of equations (pTTj) and (p3"|) . equation ([74]) simplifies to equation (|65b|) : 

^ = -2(u 2 + Q 2 ) /\{£,P}) - ^(t)(p 2 ) + 2D(t) . 

We thus conclude that we can derive the master equation for the correlators from the Kadanoff-Baym equations 
using the perturbative approximation in equations (|70[) and (|71|) . Clearly, this approximation invalidates the intricate 
rcsummation techniques of the quantum field theoretical 2PI scheme. In the 2PI framework, one resums an infinite 
number of Feynman diagrams in order to obtain a stable and thermalised late time evolution. By approximating the 
memory integrals in the Kadanoff-Baym equations, the master equation spoils this beautiful property. 

The derivation presented here can be generalised to quantum field theory. By using similar approximations, one 
can thus derive the renormalised correlator equations that would follow from the perturbative master equation. 

V. RESULTS: ENTROPY GENERATION IN QUANTUM FIELD THEORY 

Let us now return to field theory and solve for the statistical propagator and hence fix the Gaussian von Neumann 
entropy of our system. For completeness, let us here just once more recall equation (f2l)]) and (|2"5)l for the causal and 
statistical propagator: 

(d 2 + k 2 +m$)A£(M,i') - {% + k 2 )J i di 1 ^(M,fi)AS(M 1 ,f) dtiM£ th (M,*i)A£(Mi,0 = 0(76a) 

rto r* 



(% + k? + ml)F+(k,t,t') - (d 2 + k 2 ) 



[° dt 1 z;(k,t,t 1 )F l ; cc (k,t 1 ,t')+ f dtiZ^{k,t,h)F^k,ti,-e) (76b) 

J — oo J to 

f° dtiZ£(M > t 1 )AJ free (M 1 ,t')- f dhZ^k^A^tut') 

J —oo Jtn 



" dhMl th (k,t,ti)F$ ee (k, t u t')- I dtiJW£ th (M,*i)i ; *(Mi>*') 



f° dt 1 Ml th (k,t,h)Af Iee (k,t 1 ,t')+ f dhMl th (k,t,h)Al(k,h,t')=0. 

J —oo J to 



We use all self-masses calculated previously: we need the vacuum self- masses in equation (|25[) , one of the two following 
infinite past memory kernels in equation (|30[) or pip depending on the initial conditions chosen, the thermal causal 
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self-mass in (|36|) . the vacuum-thermal contribution to the statistical self- mass in equation (|43l) and finally the high 
temperature or low temperature contribution to the thermal-thermal statistical self- mass in equation (|50|) or (|48l) . 
We are primarily interested in two cases, a constant mass for our system field and a changing one: 

m<t>(t) = mo = const (77a) 
m%t) = A + Btanh(p{t-t m }), (77b) 

where we let A and B take different values. Also, t m is the time at which the mass changes, which we take to be 
pt m = 30. Let us outline our numerical approach. In the code, we take to = and we let pt and pt' run between and 
100 for example. As in the vacuum case, we first need to determine the causal propagator, as it enters the equation of 
motion of the statistical propagator. The boundary conditions for determining the causal propagator are as follows: 

A£(M) = (78a) 
d t Al(t,t%= t , = -1, (78b) 

Condition (|78a|) has to be satisfied by definition and condition (|78b|) follows from the commutation relations. 

Once we have solved for the causal propagator, we can consider evaluating the statistical propagator. As in the 
T = case, the generated entropy is a constant which can be appreciated from a rather simple argument [l|. When 
?ri0 = const, we have F,p(k, t, t') = F c j,(k, t — t') such that quantities like: 

F*(k,0) = / — W) (79a) 

POO 17,0 

Wfc,At)| At=0 = -i J — fc°W) (79b) 
d t ,d t Ft(k,At)\ At=0 = / — fc o 2 F (^), (79c) 

J -oo Z7r 

are time independent. Consequently, the phase space area is constant, and so is the generated entropy. If our 
initial conditions differ from these values, we expect to observe some transient dependence. This entropy is thus 
the interacting thermal entropy. The total amount of generated entropy measures the total amount of decoherence 
that has occurred. Given a temperature T, the thermal entropy provides a good estimate of the maximal amount of 
entropy that can be generated (perfect decoherence), however depending on the particular parameters in the theory 
this maximal amount of entropy need not always be reached (imperfect decoherence). Effectively, the interaction 
opens up phase space for the system field implying that less information about the system field is accessible to us and 
hence we observe an increase in entropy. In order to evaluate the integrals above, we need the statistical propagator 
in Fourier space: 



l iM+-{W) + iM-+{W) 
^ ' 2tM+-(k»)-iM-+(k' i ) 



kyk^ + m 2 , + iM r (k^) k^ + m 2 , + iM a (k^) 



(80) 



Here, iM r and iM a are the retarded and advanced self-masses, respectively. All the self-masses in Fourier space in this 
expression are derived in appendix [B] The discussion above is important for understanding how to impose boundary 
conditions for the statistical propagator at to. We impose either so-called "pure state initial conditions" or "mixed 
state initial conditions" . If we constrain the statistical propagator to occupy the minimal allowed phase space area 
initially, we impose pure state initial conditions and set: 

F*fa,to) = ^r— (81a) 

ZWi n 

dtF^{t,t )\ t=tQ = (81b) 
9 t ,9 t F4t,*')| t =t'= t0 = (81c) 

where uj- ln refers to the initial mass m^ito) of the field if the mass changes throughout the evolution. This yields 
Afc(io) = 1 such that: 

Sfc(to) = 0, (82) 



Initially, we thus force the field to occupy the minimal area in phase space. Clearly, if we constrain our field to be 
in such an out-of-equilibrium state initially, we should definitely not include all memory kernels pretending that our 
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field has already been interacting from negative infinity to to- Otherwise, our field would have thcrmaliscd long before 
to and could have never began the evolution in its vacuum state. If we thus impose pure state initial conditions, we 
must drop the " thermal memory kernels" : 

° dtiM£ th (M,ti)i$ w, (Mi,i') and / ° diiM^ th (fc, t, fi)A^' free (fe, t±, t') , (83) 

-oo J — oo 

but rather keep the "vacuum memory kernels" in equation (|76b[) , which are the other two memory kernels involving 
free propagators. We evaluated the relevant integrals in closed form in equation (|31[) . This setup roughly corresponds 
to switching on the coupling h adiabatically slowly at times before to- At £o ? the temperature of the environment 
is suddenly switched on such that the system responds to this change from to onwards. Note that if we would not 
include any memory effects and switch on the coupling h non-adiabatically at to, the pure state initial conditions 
would correspond to the physically natural choice. This would however also instantaneously change the vacuum of our 
theory, and we would thus need to renormalise our theory both before and after to separately. Including the vacuum 
memory kernels is thus essential, as it ensures that our evolution is completely finite at all times without the need for 
time dependent counterterms 2 . 

Secondly, we can impose mixed state boundary conditions, where we use the numerical values for the statistical 
propagator and its derivatives calculated from equations (|79")) and ([50)1 . such that we have Ak(to) = A ms = const and: 

S k (t o ) = S ias >0, (84) 

where we use the subscript "ms" to denote "mixed state" . In other words, we constrain our system initially to be in 
the interacting thermal state and S ms is the value of the interacting thermal entropy. The integrals in equation (|79[) 
can now be evaluated numerically to yield the appropriate initial conditions. For example when f3p = 1/2, k/p = 1, 
fn^l p = 1 and h/ p = 3, we find: 

F 4> (k/p=l,At)\ At=0 = 1.89885 (85a) 
d t F 4> (k/p=l,At)\ At=0 = (85b) 
d t ,d t F !t) (k/p=l,At)\ At=0 = 2.08941. (85c) 

Clearly, equation (|85b[) always vanishes as the integrand is an odd function of k°. The numerical value of the phase 
space area in this case follows from equations ([55)1 and §5§ as: 



The interacting thermal entropy hence reads: 



A ms = 3.98371 . (86) 



5^ = 1.67836. (87) 



The mixed state initial condition basically assumes that our system field has already equilibrated before to such that 
the entropy has settled to the constant mixed state value. In this case, we include of course both the vacuum memory 
kernels and the thermal memory kernels 3 . 

A few more words on the memory kernels for the mixed state boundary conditions are in order. For the vacuum 
memory kernels, we of course use equation (|30p . It is unfortunately not possible to evaluate the thermal memory 
kernels in closed form too. The two integrals in equation ([53"|) have to be evaluated numerically as a consequence. 
One can numerically verify that the integrands are highly oscillatory and do not settle quickly to some constant value 
for each t and t' due to the competing frequencies lj and k. We chose to integrate from -300 to top = and smooth 
out the remaining oscillations of the integral by defining a suitable average over half of the period of the oscillations. 

Finally, let us outline the numerical implementation of the Kadanoff-Baym equations ([76]) . Solving for the causal 
propagator is straightforward as equation ((75)) provides us for each t' with two initial conditions at t = t' and at 



2 In [70l , [7lll the renormalisation of fermions in an expanding Universe is investigated where a similar singularity at the initial time to is 
encountered. It could in their case however be removed by a suitably chosen Bogoliubov transformation. 

3 Let us make an interesting theoretical observation that to our knowledge would apply for any interacting system in quantum field theory. 
Suppose our coupling h would be time independent. Suppose also that the system field <f> and the environment field \ form a closed 
system together. Now, imagine that we are interested in the time evolution of the entropy at some finite time to- Our system field has 
then already been interacting with the environment at times before to such that one can expect that our system has equilibrated at to. 
Hence, to allow out-of-cquilibrium initial conditions, one must always change the theory slightly. The possibility that we advocate is to 
drop those memory kernels that do not match the chosen initial condition. In this way, the evolution history of our field is consistent. 
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t = t' + At, where At is the numerical step size. We can thus solve the causal propagator as a function of t for each 
fixed t' . Solving for the statistical propagator is somewhat more subtle. The initial conditions, e.g. in equation (|81l) . 
for a given choice of parameters only fix F<f,(t , to), F^(to + At, t ) = ^(to, to + At) and F^(to + At, t + At). This is 
sufficient to solve for F${t, to) and F^(t, to + At) as functions of time for fixed t' = to and t' = to + At. Now, we can 
use the symmetry relation F^(t,t') = F$(t! ,t) such that we can also find F^(to,t') and F^to + At, t') as functions of 
t' for fixed t = to and t = to + At. The latter step provides us with the initial data that is sufficient to find F,p(t,t') 
as a function of t for each fixed t'. 

Once we have solved for the statistical propagator, our life becomes much easier as we can immediately find the 
phase space area via relation ([5]). The phase space area fixes the entropy. 



A. Evolution of the Entropy: Constant Mass 

Let us firstly turn our attention to figure [5] This plot shows the phase space area as a function of time at a fairly 
low temperature ftp = 2. Starting at Afe(to) = 1, its evolution settles precisely to A ms , indicated by the dashed black 
line, as one would expect. From the evolution of the phase space area, one readily finds the evolution of the entropy 
as a function of time in figure [6] 

At a higher temperature, ftp = 1/2, we observe in figures [7] and [3] that the generated phase space area and entropy 
as a function of time is larger. This can easily be understood by realising that the thermal value of the entropy, set 
by the environment, provides us with a good estimate of the maximal amount of decoherence that our system can 
experience. Again we observe an excellent agreement between A ms or S^s and the corresponding numerical evolution. 

Let us now discuss figure [9l Here, we show two separate cases for the evolution of the entropy: one at a very low 
temperature ftp = 10 (in black) and one vacuum evolution ftp = oo (in gray) which we already calculated in [lj. As 
we would intuitively expect, we see that the former case settles to an entropy S ms = 0.04551 that is slightly above 
the vacuum asymptote S ms = 0.04326. 

Finally, in figure [TU] we show the interacting phase space area A ms as a function of the coupling h. For h/p <C 1, 
we see that A ms approaches the free thermal phase space area Af rec = coth(/3w/2). For larger values of the coupling, 
we see that A ms > Af reo . If these two differ significantly, we enter the non-perturbative regime. In the perturbative 
regime, this plot substantiates our earlier statement that the free thermal entropy Af ree provides us with a good 
estimate of the total amount of decoherence that our system can experience. Our system however thermalises to A ms , 
and not to Af ree as the interaction changes the nature of the free thermal state. 

The most important point of the results shown here is that, although a pure state with vanishing entropy Sk = 
remains pure under unitary evolution, we perceive this state over time as a mixed state with positive entropy S ms > 
as non-Gaussianities are generated by the evolution (both in the correlation between the system and environment as 
well as higher order correlations in the system itself) and subsequently neglected in our definition of the Gaussian von 
Neumann entropy. The total amount of decoherence corresponds to the interacting thermal entropy S ma . 



B. Decoherence Rates 



As the Gaussian von Neumann entropy in equation (|4|) is the only invariant measure of the entropy of a Gaussian 
state, we take the point of view that this quantity, or equivalently the phase space area in equation ([5]), should be 
taken as the quantitative measure for decoherence. This agrees with the general view on decoherence according to 
which the decoherence rate is the rate at which a system in a pure state evolves into a mixed state due to its interaction 
with an environment. This is to be contrasted with some of the literature where different, non- invariant measures are 
proposed (fj [l5| . For example in [TBI . [72j , the superposition of two minimum uncertainty Gaussian states located at 
positions x and x' is considered. The decoherence rate is defined differently, i.e.: it is the characteristic timescale at 
which the off-diagonal contributions in the total density matrix decay and coincides with the timescale at which the 
interference pattern in the Wigner function decays. It is given by: 



r^ 1 



At 



where the thermal de Broglie wavelength is given by At = (2m/csT)~ 1 / 2 . In other words, according to 15|, the 
decoherence rate depends on the spatial separation a; — a;' of the two Gaussians. Note that in quantum field theory the 
expression would generalise to 1 cx (0 — (j)') 2 - This is just one example, one can find other definitions of decoherence 
in the literature. 
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Figure 5: Phase space area as a function of time. It settles 
nicely to A ms , indicated by the dashed black line. We use 
f3p = 2, k/p = 1, m^/ p — 1, h/p = 4 and a totaZ number 
of steps N = 2000 up to tp = 100. 
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Figure 6: Entropy as a function of time. The evolution of 
the entropy is obtained from the phase space area in figure 
[3[ For 50 < tp < 100, the entropy continues to coincide 

With Sms ■ 
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Figure 7: Phase space area as a function of time. At 
high temperatures, we see that the phase space area set- 
tles quickly again to A ms . We use ftp = 1/2, k/p = 1, 
m^/p = 1, h/p = 3 and N — 2000 up to tp = 100. 



Figure 8: Entropy as a function of time, which follows 
again from the evolution of the phase space area as a func- 
tion of time as depicted in figure^ Our pure state quickly 
appears to our observer as a mixed state with a large pos- 
itive entropy Sms- 
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Figure 9: Entropy as a function of time. The evolution 
at very low temperatures ftp — 10 (black) resembles the 
vacuum evolution f3p = oo (gray) from |j/ as one would 
intuitively expect. We furthermore use k/p = 1, m^/p = 
1, h/p — 4 and N = 2000 up to tp = 100. 



Figure 10: We show A ms , the interacting thermal phase 
space area, as a function of h/p. For h/p <C 1, we see that 
A ms is almost equal to the free phase space area Af reB = 
coth(/3aj/2) indicated by the dashed black line. For larger 
values of h/p we approach the non-perturbative regime. 
We use ftp = 1/2, k/p — 1 and m^/ p — 1. 
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The main difference is that our decohcrencc rate docs not depend on the configuration space variables x or <f> but is 
an intrinsic property of the state. In other words, we do not look at different spatial regions of the state, but rather 
to the state as a whole from which we extract one decoherence rate. As we outlined in Q, a nice intuitive way to 
visualise the process of decoherence is in Wigner space. The Wigner transform of a density matrix coincides with 
the Fourier transform with respect to its off-diagonal entries. As discussed previously, the phase space area measures 
the area the state occupies in Wigner space in units of the minimum phase space area ft/2, which we refer to as the 
statistical particle number n. The pure state considered in the previous subsection decoheres and its phase space area 
increases to approximately its thermal value. When A ^> 1 (n » 1), different regions in phase space of area ft/2 are, 
to a good approximation, not correlated and thus evolve independently. As we have considered Gaussian states only 
and not the superposition of two spatially separated Gaussians, which when considered together is in fact a highly 
non-Gaussian state, a direct comparison is not straightforward. 

Let us extract the decoherence rate from the evolution of the entropy. We define the decoherence time scale to be 
the characteristic time it takes for the phase space area Afc(t) to settle to its constant mixed state value A ms . The 
phase space area approaches the constant asymptotic value in an exponential manner: 

^5A k (t)+T dec 5A k (t)=0, (89) 

where 5Ak(t) = A ms — Afc(t) and where Tdoc is the decoherence rate. This equation is equivalent to hk = —Tdcc(nk — 
n ms ), where is defined in equation (J5J and n ms is the stationary n corresponding to A ms . As in the vacuum case 
we anticipate that the decoherence rate is given by the single particle decay rate of the interaction <f> —> \ 2 ■ The 
single particle decay rate reads 4 : 



Im(.MJ) 



W0 



lo g I ; - „, — rr . (9°) 



where we used the retarded self-mass in Fourier space in equation (|B3ap and several relevant self-masses in appendix |B| 
Let us briefly outline the steps needed to derive the result above. In order to calculate zM^(fc M ), we use tMf^.Qt 1 *) 

in equation (jB6a| and iMf^ c (k ,i ) in equation (|B12|1 . There are no thermal-thermal contributions to iM^(fc^) which 
can be appreciated from equation (|B10[) . Finally, in order to derive the vacuum-thermal contribution, let us recall 
equation (|24c|) given by: M+ + (k,t,t') = M£(k,t,t') + sgn(t-t')iM%(k,t,t')/2. We clearly need the vacuum-thermal 
contribution to Mf which is given in equation (|B 14|) . The imaginary part of the second term vanishes, which can 
be seen by making use of an inverse Fourier transform, just as in the first lines of equations (|B17|) and (|B18|) . This 
fixes «M^(fc M ) completely. 

One should calculate the imaginary part of the retarded self-mass as it characterises our decay process, which 
follows from equation ([80]). In order to calculate the decay rate, we have to project the retarded self- mass on the 
quasi particle shell k° = oja. Of course, one should really take the perturbative correction to the dispersion relation 
of order C(ft 2 /w^) into account but this effect is rather small. Alternatively, we can project the advanced self- mass 
in Fourier space on fc° = — oj$. We thus expect: 



dec 



4>^xx ■ (91) 



Let us examine figures [TT1 and [T2l From our numerical calculation, wc can thus easily find 6Ak(t) = A ms — Ak{t) 
which we show in solid black. We can now compare with the single particle decay rate in equation (|90[) and plot 
— Tj,-> xx t. We conclude that the decoherence rate can be well described by the single particle decay rate in our model, 
thus confirming equation (|91|) above. 



C. The Emerging fe° = Shell 



To develop some intuition, we depict F^,(k fM ) as a function of k° keeping various other parameters fixed. In the 
vacuum j3p — oo, it is clear from the analytic form of the statistical propagator that a k° = shell does not exist. In 
the vacuum, we have that F^ > (k fJl ) = for < k. At low temperatures, /3p = 2, we observe in figurcfT3lthat two more 
quasi particle peaks emerge where < k. The original quasi particle peaks at |fc°| ~ however still dominate. 



4 For cases where m x ^ 0, see [73ll . 
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Figure 11: Decoherence rate at low temperatures. We 
show the exponential approach to A ms in solid black and 
the corresponding decoherence rate given in equation 
$91\) (dashed line). We use the phase space area from 
figure [3[ 
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Figure 12: Decoherence rate at high temperatures. We 
show the exponential approach to A ms in solid black and 
the corresponding decoherence rate given in equation 
S91]) (dashed line). We use the phase space area from 
figure^ 



F(k") 



m 




60 




MO- 

M J 
l 1 // 


V " 

V i \ 


u> J 20 


V 1 1 
A 1 ' « 




/ \y \ 

~ " ' 


~*r^ / ^ 





-3-2-10 1 2 3 



Figure 13: The statistical propagator at low temperatures. 
Apart from the original quasi particle peaks, two more 
peaks emerge at |fe°| < k. We use f3p = 2, k/p — 1, 
iritjy/p = 1 and h/p — 4 (solid black), h/p = 2 (dashed). 
In the latter case, we do not show the entire original quasi 
particle peak for illustrative reasons. 



Figure 14: The statistical propagator at high tempera- 
tures. For larger and larger coupling, we observe that 
the newly emerging quasi particle peaks move closer to 
k° = 0, where they eventually almost completely over- 
lap. We use ftp = 0.1, k/p — 1, m^/p = 1 and 
h/p — 1.61 (solid black), h/p — 1.5 (dot-dashed) and 
h/p = 1 (dashed). 



At high temperatures, j3p = 0.1, we observe in figure [Ml that the two additional quasi particle peaks already present 
at lower temperatures increase in size and move closer to k° = 0, where they overlap. The original quasi particle 
peaks located at |fc°| ~ broaden as the interaction strength h increases. Moreover, for increasing h, the original 
quasi particle peaks get dwarfed by the new quasi particle peaks at | A: 1 < k that by now almost completely overlap 
at k° = 0. 

What we observe here is related to the coherence shell at k° — first introduced by Herranen, Kainulainen and 
Rahkila [13, 111] to study quantum mechanical reflection and quantum particle creation in a thermal field theoretical 
setting (and for a discussion of fermions see [32|, HH). They interpret this new spectral solution of the statistical two 
point function as a manifestation of non-local quantum coherence. As we have just seen, the statistical propagator 
at late times will basically evolve to equation (|5U]) . We conclude that the emerging fc° = shell translates to large 
entropy generation at high temperatures. It is also clear that the naive quasi particle picture of free thermal states 
breaks down in the high temperature regime. 
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D. Evolution of the Entropy: Changing Mass 



Let us now study the evolution of the entropy where the mass of the system field changes according to equation 
(|77b[) . For a constant mass ms, the statistical propagator depends only on the time difference of its arguments 
Fj,(k, t, t') = F c f,(k, t — t') due to time translation invariancc. This observation allowed us to find the asymptotic value 
of the phase space area by means of another Fourier transformation with respect to t — t' . When the mass of the 
system field changes, however, we introduce a genuine time dependence in the problem and we can only asymptotically 
compare the entropy to the stationary values well before and after the mass change. It is important to appreciate 
that the counterterms introduced to renormalise the theory do not depend on so we do not have to consider 
renormalisation again 

Depending on the size of the mass change, we can identify the following two regimes: 

\Pk\ ^ 1 adiabatic regime (92a) 

|/9/c| 2 ^ 1 non — adiabatic regime, (92b) 

where ftk is one of the coefficients of the Bogoliubov transformation that relates the initial (in) vacuum to the final 
(out) vacuum state. As a consequence of the mass change, the state gets squeezed Q. If ftu = 0, the in and out 
vacuum state are equal such that \ftk\ 2 quantifies the amount of particle creation and reads [75| : 

sinh 2 ( ^«p 



sinh (2^) sinh f 2nuj ™ V 2 P J 

+ k 2 are the initial and final frequencies. Also, m? in = A — B and 
TO | out = ^ + 1? where we made use of equation ([77bJ) . Finally, we defined uj± = |(w ou t ±w. n ). The word "particle" in 
particle creation is not to be confused with the statistical particle number defined by means of the phase space area in 
equation ([§]) . Whereas the latter counts the phase space occupied by a state in units of the minimal uncertainty wave 
packet, the former corresponds to the conventional notion of particles in curved spacetimes where one plane wave 
field excitation at|0) = \k) is referred to as one particle (for a discussion on wave packets in quantum field theory, 

see [zl,[z3). When we consider a changing mass in the absence of any interaction terms, \ftk\ 2 increases whereas the 
phase space area remains constant. For the parameters we consider in this paper \ft k \ 2 ~ O(10 -4 ) such that we are 
in the adiabatic regime. 

Let us consider the coherence effects due to a mass increase and decrease in figures [T__ HH EZl and [T_] Here, we 
take m^/p = 1 and m^/p = 2 giving rise to the constant interacting thermal entropies 5ms and Sms, respectively. 
The numerical value of these asymptotic entropies is calculated just as in the constant mass case such that we find 
Sms > Sml- We use mixed state initial conditions as outlined in equation and moreover we insert the initial 
mass in the memory kernels. In figure [15] we show the effects on the entropy for a mass increase at fairly low 
temperatures ftp = 2. In gray we depict the two corresponding constant mass entropy functions to compare the 
asymptotic behaviour. In order to calculate the latter, we also use mixed state boundary conditions. Clearly, well 

before and after the mass increase, the entropy is equal to the constant interacting thermal entropy, Sm£ and Sml, 

(2) 

respectively. The small difference between the numerical value of the interacting thermal entropy S m i (in dashed 
gray) and the corresponding m^/p = 2 constant mass evolution is just due to numerical accuracy. It is interesting to 
observe that the new interacting thermal entropy is reached on a different time scale than the one at which the 
system's mass has changed. Again, we verify that the rate at which the phase space area changes, defined analogously 
to equation ([89]) . can be well described by the single particle decay rate ([90]). Given the fact that the mass changes so 
rapidly in our case, one should use the final mass m^. out in equation (|90p . In figure [T9l we show both the exponential 

( 2) _____ 

approach towards the constant interacting phase space area A m g and the decay rate (|9T)]) . In order to produce figure 
1191 we subtract the constant mass evolution of the phase space area using mixed state initial conditions rather than 
A^l to find SA k (t) in equation ([89]) . 

This qualitative picture does not change when we consider the same mass increase only now at higher temperatures 
ftp = 1/2 in figure [TBI The interacting thermal entropies in this case are larger due to the fact that the temperature is 

(2) 

higher. Again we observe a small difference between S m £ and the m^/p = 2 constant mass evolution due to numerical 
accuracy. Also, the decoherence rate can be well described by the single particle decay rate which we depict in figures 
(__]) and (2_J. 

When we consider the "time reversed process", i.e.: a mass decrease from m^/p = 2 to m^/p = 1, we observe 
an entropy increase. We show the resulting evolution of the entropy in figures [T71 and [T51 for ftp = 2 and ftp = 1/2, 
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Figure 15: Entropy as a function of time for a mass in- 
crease from m^/p = 1 to m^/p = 2, giving rise to the 
constant interacting thermal entropies si! s ' and Smi, re- 
spectively. The mass changes rapidly at tp = 30. We use 
I3p = 2, k/p = 1, h/p — 4 and N — 1600. The gray 
lines are the corresponding constant mass entropy func- 
tions where we use mixed state initial conditions. 
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Figure 16: Entropy as a function of time for a mass in- 
crease from m^/p — 1 to m^/p = 2, giving rise to the 
constant interacting thermal entropies Sm2 and sS, re- 
spectively. The mass changes rapidly at tp = 30. We use 
Pp = 1/2, k/p = 1, h/p — 3 and N = 1600. The gray 
lines are the corresponding constant mass entropy func- 
tions where we used mixed state initial conditions. 
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Figure 1 7: Entropy as a function of time for a mass de- 
crease where we used the same parameters as in fiaure [T5l 



Figure 18: Entropy as a function of time for a mass de- 
crease where we used the same parameters as in Haure \16\ 



respectively. The evolution of the entropy reveals no further surprises and corresponds to the time reversed picture of 
figures [TS] and [Tf)J The decoherence rate for a mass decrease can again be well described by the single particle decay 
rate in equation (|90|) . 

We observe that the rate at which the mass changes is much larger than the decoherence rate. As long as this 
condition is satisfied, coherence effects continue to be important. Eventually though, the Gaussian von Neumann 
entropy settles to its new constant value and no particle creation remains as our state thermalises again. In the 
context of baryogenesis, we thus expect that quantum coherence effects remain important as long as this condition 
persists too. Of course, one would have to generalise our model to a CP violating model in which the effects that are 
of relevance for coherent baryogenesis scenarios are captured. 



E. Squeezed States 



The effect of a large non-adiabatic mass change on the quantum state is a rapid squeezing of the state which can 
neatly be visualised in Wigner space. Although it is numerically challenging to implement a case where the mass 
changes non-adiabatically fast, we can probe its most important effect on the state by considering a state that is 
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Figure 19: Decoherence rate at low temperatures. We 
show the exponential approach to A„] in solid black and 
the corresponding decoherence rate given in equation f91p 
line). We use the phase space area from fiaure[To[ 



Figure 20: Decoherence rate at high temperatures. We 
show the exponential approach to in solid black and 
the corresponding decoherence rate given in equation H91\) 
(dashed line). We use the phase space area from fiaure lTSl 



significantly squeezed initially. A pure and squeezed state is characterised by the following initial conditions: 

F c f,(k, t , to) = - — [cosh(2r) — sinh(2r) cos(2<y9)] 



d t drF4k,t,t% =t , =to 
d t F,p{k,t,to)\ 



2w v 



cosh(2r) + sinh(2r) cos(2<^)] 



sinh(2r) sm(2ip) . 



(94a) 
(94b) 
(94c) 



Here, <p characterises the angle along which the state is squeezed and r indicates the amount of squeezing. As a 
squeezed state is pure, we have A&(io) = 1 initially. A mixed initial squeezed state condition can be achieved by 
multiplying equation (|94|) by a factor. 

We show the corresponding evolution for the phase space area in two cases in figures [5T] and [521 As the squeezed state 
thermalises, we observe two effects. Firstly, there is the usual exponential approach towards the thermal interacting 
value A ms we observed before. As we showed previously, this process is characterised by the single particle decay rate 
in equation (|90[) . Secondly, superimposed to that behaviour, we observe damped oscillatory behaviour of the phase 
space area as a function of time that is induced by the initial squeezing. 

The latter process in principle introduces a second decay rate in the evolution: one can associate a characteristic 
time scale at which the amplitude of the oscillations decay (superimposed on the exponential approach towards A ms ). 
One can read off from figures [5T] and [521 that the exponential decay of the envelope of the oscillations can also be 
well described by the single particle decay rate in equation (fUU|) . We thus observe only one relevant time scale of the 
process of decoherence in our scalar field model: the single particle decay rate. We thus conclude that in the case of 
a non-adiabatic mass change, the decay of the amplitude of the resulting oscillations will be in agreement with the 
single particle decay rate too. 



VI. CONCLUSION 



We study the decoherence of a quantum field theoretical system in a rcnormaliscd and perturbativc 2PI scheme. As 
most of the non-Gaussian information about a system is experimentally hard to access, we argue in our "correlator 
approach" to decoherence that neglecting this information and, consequently, keeping only the information stored in 
Gaussian correlators, leads to an increase of the Gaussian von Neumann entropy of the system. We argue that the 
Gaussian von Neumann entropy should be used as the quantitative measure for decoherence. 

The most important result in this paper is shown in figure [5J where we depict the time evolution of the Gaussian 
von Neumann entropy for a pure state at a high temperature. Although a pure state with vanishing entropy Sk = 
remains pure under unitary evolution, the observer perceives this state over time as a mixed state with positive 
entropy 5 ms > 0. The reason is that non-Gaussianitics arc generated by the unitary evolution (both in the correlation 
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Figure 21: Phase space area as a function of time for a 
squeezed initial state. We use tp = 0, e 2r = 1/5, f3p = 0.5, 
k/p — 1, m^/p — 1, = 3 and iV = 300 tip to tp = 15. 
27ie gray line indicates the pure state evolution previously 
considered in figure [7] 



Figure 22: Phase space area as a function of time for 
a squeezed initial state. On top of the usual exponential 
approach towards A ms , we observe oscillatory behaviour. 
The amplitude of the oscillations decays with the single 
particle decay rate as well. We use tp = 0, e 2r = 5 and 
the other parameters are given in figure rjJJ 



between the system and environment as well as in higher order correlations in the system itself) and subsequently 
neglected in our Gaussian von Neumann entropy. 

We have extracted two relevant quantitative measures of decoherence: the maximal amount of decoherence S ms 
and the decoherence rate Tdoc- The total amount of decoherence corresponds to the interacting thermal entropy S ms 
and is slightly larger than the free thermal entropy, depending on the strength of the interaction h. The decoherence 
rate can be well described by the single particle decay rate of our interaction T^ xx . 

This study builds the quantum field theoretical framework for other decoherence studies in various relevant situations 
where different types of fields and interactions can be involved. In cosmology for example, the decoherence of scalar 
gravitational perturbations can be induced by e.g. fluctuating tensor modes (gravitons) p, isocurvature modes [go[ 
or even gauge fields. In quantum information physics it is very likely that future quantum computers will involve 
coherent light beams that interact with other parts of the quantum computer as well as with an environment [78|, [T9| . 
For a complete understanding of decoherence in such complex systems it is clear that a quantum field theoretical 
framework such as developed here is necessary. 

We also studied the effects on the Gaussian von Neumann entropy of a changing mass. The Gaussian von Neumann 
entropy changes to the new interacting thermal entropy after the mass change on a time scale that is again well 
described by the single particle decay rate in our model. It is the same decay rate that describes the decay of 
the amplitude of the oscillations for a squeezed initial state. One can view our model as a toy model relevant for 
electroweak baryogenesis scenarios. It is thus interesting to observe that the coherence time scale (the time scale at 
which the entropy changes) is much larger than the time scale p _1 at which the mass of the system field changes. We 
conclude that the coherent effect of a non-adiabatic mass change (squeezing) does not get immediately destroyed by 
the process of decoherence and thermalisation. 

Finally, we compared our correlator approach to decoherence to the conventional approach relying on the perturba- 
tive master equation. It is unsatisfactory that the reduced density matrix evolves non-unitarily while the underlying 
quantum theory is unitary. We are not against non-unitary equations or approximations in principle, however, one 
should make sure that the essential physical features of the system one is describing are kept. The perturbative 
master equations does not break unitarity correctly, as we have shown in this paper. On the practical side, the master 
equation is so complex that field theoretical questions have barely been addressed: there does not exist a treatment 
to take perturbative interactions properly into account, nor has any reduced density matrix ever been renormalised. 
This is the reason for our quantum mechanical comparison, rather than a proper field theoretical study of the re- 
duced density matrix. In section llVDI however, we outline the perturbative approximations used to derive the master 
equation from the Kadanoff-Baym equations, i.e.: in the memory kernels of the Kadanoff-Baym equations we insert 
free propagators with appropriate initial conditions. A proper generalisation to derive the renormalised perturbative 
master equation in quantum field theory from the Kadanoff-Baym equations should be straightforward. In the simple 
quantum mechanical situation, we show that the entropy following from the perturbative master equation generically 
suffers from physically unacceptable secular growth at late times in the resonant regime. This leads to an incorrect 
prediction of the total amount of decoherence that has occurred. We show that the time evolution of the Gaussian 
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von Neumann entropy behaves well in both the resonant and in the non-resonant regime. 
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Appendix A: Derivation of M^ th _ th (k, At) 



Only the high and low temperature limits of M^ th _ th (k, At) can be evaluated in closed form. We derive these 
expressions in this appendix. 



1. Low Temperature Contribution 



Let us recall equation (|39b[) where we can perform the w-intcgral by making use of equation (|3"3")l and l/(e^ w — 1) 

Eoo 
n=l' 



Mf th _ th (A,At) 



h 2 r „ cos(fc 1 At) ^ 
1 e ^fci _ i 



llT 2 k 



-n(3uj 







^ (At) 2 + (npf 



-n(3cos(ujAt) + Aisin(wAt)] , (Al) 



where ui\ — (k ± ki) 2 + m 2 . We now prepare this expression for hi integration by making use of l/(e° 1 — 1) 
Sm=i e _/3mfcl and some familiar trigonometric identities: 



t(At) 2 + W) 3 io 1 



16n 2 k ^ (At) 

m,n—l 



(A2) 



x < — (3ne 



-f3n(k+ki) 



cos[(2fci + fc)At] + cos(fcAt) 



sin[(2fci + A;)Ai] + sin(fcAt) 



+ Ate~ /3,l(fc+fcl) 

+6(k-k 1 )(3ne- l3ni - k - kl) cos[(2fti-jfc)At] + cos(fcAt)j + 9(k-k 1 )At e -M*-*i) sin[(2ifci - k)At] - sin(fcAi) 
+6»(fci - fc)^n e -^( fc i- fc ) [ cos[(2fci - fc) At] + cos(fcAt)] - 9(ki - k)At e -£"Oi-*0 \ s i n [(2A; 1 - fc) Ai] - sin(fcAt) 



Upon integrating over k\ and rearranging the terms we obtain: 

1 



167r 2 fc ^ (At) 2 + (n/3) 5 
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sin(fcAi) ( 
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At 



At 



/3At(m + 3n) 
m + n)] 2 + (2Ai) 2 j3(m-n) 

(A3) 



3(fcAi)(e 



[/3(m - n)] 2 + (2Ai) 2 /?(m + n) 

-/3 2 n(n + m) + 2(At) 2 /3 2 n(m - n) + 2(At) 2 n n 
[f3(m + n)] 2 + (2At) 2 [j3(m - n)} 2 + (2At) 2 m + n m-n 



This expression contains two singular terms when m = n. By performing the integral (|A2[) in that case, they are to 
be interpreted as: 



-Bnk ~ — Bmk 

e ra—n 



m — n 



13k e 



-fink 



(A4) 



This expression allows us to obtain the low temperature /3k 3> 1 limit of M^ th _ th (k, At). It then suffices to consider 
three contributions in equation (|A3[) only. Firstly, there is the contribution for m = 1 = n, for n = 1 and m > 2, and 
finally for m = 1 and n > 2. The sum in the last two cases can be evaluated in closed form, such that one obtains 
equation (|48l) . 
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2. High Temperature Contribution 

Let us now consider the high temperature limit. It is clear from equation (|A3|) that when /3A; <SC 1 there is 
unfortunately no small quantity to expand about as both m and n can become arbitrarily large. Therefore, we 
go back to the original expression (|39b[) . proceed as usual by making use of (|33p and rewrite it in terms of new 
(u, u)-coordinates ("lightcone coordinates"), defined by: 

u = k\ — lu (A5a) 
v = ki+u, (A5b) 

such that of course k\ = (v + u)/2 and uj = (v — u)/2, in terms of which the region of integration becomes: 

- k < u<k (A6a) 
k < v < oo. (A6b) 

Equation (|39b[) thus transforms into: 



Mi 
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(k,At) 



du 



32ir 2 k J_ k 2sinh(/3w/2) J k 



dv [cos(uAt) + cos(vAt)] 



(v+u)/2 _ I c P(v-u)/2 _ I 



(A7) 



where we took account of the Jacobian J = \d(ki,u>)/d(u,v)\ = 1/2. One can now perform the w-integral involving 
the cos(uAi)-term. Secondly, since we are interested in the limit /3k -C 1 and we moreover have \u\ < k, note that we 
also have \/3u\ <C 1. The cos(tAi)-term can thus be expanded around |/3u| <C 1. An intermediate result reads: 
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(A8) 



The reader can easily see that we deliberately not Taylor expand the second line fully around \/3u\ <C 1. The reason 
is that the subsequent integration renders such a naive Taylor expansion invalid. Let us first integrate the first line 
of equation (|A8[) . We now expand the cos(itAi)-integral around \(3(k ±«)| < 1: 



du cos(wAi) 



log (1 - exp[-/?(fc + u)/2]) log (1 - exp[-/3(fc - u)/2\) 
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Let us first evaluate the simple it-integrals in equation (|A9[) : 
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where ci(z) and si(z) are the cosine and sine integral functions, respectively, defined in equation 
complicated integral in (|A9[) is: 
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By making use of its definition, we expand the hypergeometric function \F 2 as follows: 
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Inserting this into (|A11[) and performing the n-sum we obtain: 
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where we performed the n-sum for m — separately. Note that: 
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Finally, we can perform the m-sum appearing in equation (|A13|) to yield: 
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(A15) 

It is useful to know the expansions of the hypergeometric function in ()A15|) . The large time (kAt 3> 1) expansion of 
this function is: 
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whereas the small times (kAt <C 1) limit yields: 
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We still need to perform some more integrals in equation 
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We can now perform the u-integral: 
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d V e-^ (n+2)+ "' At , 



where the reader can easily verify that the argument of both logarithms in the first line is positive. In the second 
line we have expanded the logarithm and made use of the binomial series. Due to the cosine appearing in equation 
(|A18[) . we are only interested in the real part of the integral on the second line. The w-integral can now trivially be 
performed. In order to extract the high temperature limit correctly, it turns out to be advantageous to perform the 
m-sum in equation (|A17|) first: 
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The Hypergeometric function can be expanded in the high temperature limit as: 
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We have checked using direct numerical integration that the analytic answers improves much if we keep also the 
second order term in this expansion. Finally, we can perform the remaining sum over n, yielding: 
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where of course we are interested in the real part of the expression above. Having performed all the integrals needed 
to calculate the high temperature limit of M^ tYl _ th {k, At), we can collect the results in equations (|A8[) . (|A10[) . (| A15I) . 
(|A17j) and equation (|A21|) above, finding precisely equation (|50[) . 



Appendix B: The Statistical Propagator in Fourier Space 



lAt-ffc") = 9 — , (Bib) 



This appendix is devoted to calculating the statistical propagator in Fourier space at finite temperature. The two 
Wightman functions, needed to calculate the statistical propagator through equation (jllbj) . are given by: 

iAT + (kn = * , * ' (Bla) 

where we have made use of the definition of the advanced propagator: 

= W+nt + tMl^m ' (B2) 

and the definitions of the advanced and retarded self-masses: 

iM; tIen (k») = iM++ n (kn - iM+- (W) = *M^ + (fc") - *M^(fc") (B3a) 

%M$ ttea (W) = d/++„(F) - iM^+(fc") = iM+-(H>) - iM^- n (k») , (B3b) 

Our starting point is: 
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The thermal propagators appearing in this equation are of course given by (|19|) . where m x — > 0. This calculation 
naturally splits again into three parts: 
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where: 
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^t + h-th(^) = ~ J ^p 47 ^ " k 'n( k , ~ K)) 6 WK) ^ q (l fc ° " fc'°IK q (l fc '°l) , (B6c) 
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where the vacuum contribution (|B6a[) has already been evaluated and renormalised in [l[ . As all thermal contributions 
are finite, we can safely let D — > 4 and make use of equation 
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Here, fc = |fc| as before. Transforming to (u, w)-coordinates already used in equation (|A5j) now yields: 
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7 2 /-fc /-oo /i \ 



fc°± -(u + v) 



[5(k° ±v) + S(k° ±u)](B8b) 



Let us firstly calculate l M^ h _ th (k fl ). The Dirac delta-functions trivially reduce equation (|B8b[) further and moreover, 
we can make use of: 
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The final result for iM^ h _ th (k^) thus reads: 
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Since this contribution to iM^ + {k^) does not depend on the pole prescription, it completely fixes similar contribu- 
tions to the other self-masses, e.g. i^tth-fiS^^) = ? ^<Ath-th(^' i )' It turns out that equation (|B8a[) is not most 
advantageous to derive iM^ &c _ th (k^) . 

Let us therefore firstly evaluate zM ±=F (fc A1 ). Let us thus start just as in equation (|B4|) and set: 

tM}*m = *M±T c (fc") + iM±^_ th (k») + iM±l_ th (kn . (Bll) 

The vacuum-vacuum contribution has been evaluated in [l[ and is given by: 

^ac^') = -^ d ^ k ° fc ) ■ (B12) 

The thermal-thermal contributions are given above in equation (jBlOj) . so we only need to determine the vacuum- 
thermal contributions. Hence, we perform an analogous calculation as for iM ++ and transform to the familiar 
lightcone coordinates u and v to find the following intermediate result: 
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The delta functions allow us to perform one of the two integrals trivially. The remaining integral can also be obtained 
straightforwardly: 
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By subtracting and adding the above self-masses, we can obtain the vacuum-thermal contributions to the causal and 
statistical self- masses in Fourier space from equations Q32[) and (|37p . respectively. The vacuum-thermal contribution 
to the causal self-mass reads: 



M£ vac _ th (P) = *M+ TOC _ th (F) - »M -+ c _ th (F) 



(B15) 
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where we have made use of the theta functions to bring this result in a particularly compact form. Likewise, the 
vacuum-thermal contribution to the statistical self-mass now reads: 
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(B16) 



167rfc/3 



sgn(fc - |fc°|) log (l - c-^ fe +l fc0 l>) + log (l - c-^l fc -l fc0 ll) 



As a check of the results above, we performed the inverse Fourier transforms of the causal and statistical self-masses 
in equations (|3"rJ]) and (gOJ), respectively, and found agreement with the results presented above. 

The most convenient way of solving the vacuum-thermal contribution to iM^ + (k^) is by making use of equation 
p4| and (|36|) . Let us set the imaginary part of iM~f + {k^) equal to: 
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where we have taken the high temperature limit on the second line. We thus have: 
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where have performed the remaining integral on the second line straightforwardly. Using equation 
find: 
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(B18) 
and (|B16|) we 

(B19a) 
(B19b) 



167rfc/3 



sgn(ft — |ft |) log ( 1 — e 2 



-|(fc+|fc°|) 



)+log( 



1 - e 



iM 



^>,vac— th 



(kn 



,3fc<i ^/l 2 



167rfc/3 



+ _ [ sgn ( fc o + k ) + sgn (-fcO + fc)] 



sgn(fc - |fc°|) log (l - c-^ fc+ l fc0 D) + log (l 



-^-[sgn(ft + fc)+sgn(-ft° + ft)] 



\k-\k 



■") 



(B20a) 



(B20b) 



One can check equation (|B18|) by means of an alternative approach. The starting point is the first line of equation 
(|3"rj)) and one can furthermore realise that differentiating Af s t+(ft' i ) with respect to fc° brings down a factor of lAt 
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which conveniently cancels the factor of At that is present in the denominator. One can then integrate the resulting 
expressions (introducing e regulators and UV cutoffs where necessary) confirming expression (|B18|) . 
In the low temperature limit, equation (|B17[) reduces to: 



M++(fc,At) ^ -^sin(fcAt)sgn(At) 



We can introduce an e regulator: 
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Analogously, we can derive the following expressions for the vacuum-thermal contributions to iA/^ :± (fc AI ) in the low 
temperature limit: 
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We have now all self-masses at our disposal necessary to calculate F(k ti ). To numerically evaluate the integrals in 
equation ([79| we do not rely on the high and low temperature expressions in equations (|B23p or (|B20j) but we rather 
use exact numerical methods, i.e.: the first line of equation (|B17p . 
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